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Abstract 

A generic Hamiltonian, which incorporates the effect of the orbital contrac- 
tion on the hopping amphtude between the nearest sites, is studied both an- 
alytically at the weak coupling limit and numerically at the intermediate and 
strong coupling regimes for finite atomic cluster. The effect of the orbital con- 
traction due to hole localization at atomic sites is specified with two coupling 
parameters V and W (multiplicative and additive contraction terms). The 
singularity of the vertex part of the two-particle Green's function determines 
the critical temperature T^ and the relaxation rate r(T) of the order param- 
eter at temperature above T^. Unlike in conventional BCS superconductors, 
r has a non-zero imaginary part which may influence the fluctuation con- 
ductivity of superconductor above Tc. We compute the ground state energy 
as a function of the particle number and magnetic flux through the cluster, 
and show the existence of the parity gap A appearing at the range of system 
parameters consistent with the appearance of Cooper instability. Numeric 
calculation of the Hubbard model (with [7 > 0) at arbitrary occupation does 
not show any sign of superconductivity in small cluster. 



I. FORMULATION OF THE MODEL 

High temperature superconductivity in the lanthanumB, yttrium! and related copper- 
oxide compounds remains a subject of intensive investigation and controversy. It was 
suggested that electron-phonon interaction mechanism, which is very successful in under- 
standing of conventional ("low temperature") superconductors within the Bardeen-Cooper- 
Schriffer scheme!, may not be adequate for high-T^ cuprates, and even the conventional 
Fermi liquid model of metallic state may require reconsideration. This opens an area for 
investigation of mechanisms of electron-electron interaction which can be relevant in under- 
standing peculiarities of superconducting, as well as normal state, properties of cuprates. 
Specific to all of them is the existence of oxide orbitals. Band calculation^'! suggest that 
hopping between the oxygen "PxiVy orbitals and between the copper d^-^^yi orbitals may be of 
comparable magnitude. On the experimental side, spectroscopic studiesOQ clearly show that 
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the oxygen band appears in the same region of oxygen concentration in which superconduc- 
tivity in cuprates is the strongest. Therefore there exists a possibihty that specific features 
of oxide compounds may be related to oxygen-oxygen hopping, or to the interaction between 
the copper and the rotational — Py collective modes. If the oxygen hopping is significant 
then it immediately follows that intrinsic oxygen carriers {px,Py oxygen holes) should be 
different from the more familiar generic s-orbital derived itinerant carriers. The difference 
is related to low atomic number of oxygen such that removing or adding of one electron to 
atom induces a substantial change in the Coulomb field near the remaining ion and there- 
fore results in the change of the effective radius of atomic orbitals near the ion. This will 
strongly influence the hopping amplitude between this atom and the atoms in its neighbor- 
hood. Such "orbital contraction" effect represents a source of strong interaction which does 
not simply reduce to the Coulomb (or phonon) repulsion (or attraction) between the charge 
carriers. It was suggested by Hirsch and coauthorsi'i0, and by the present authors0'0'E10 
that the occupation dependent hopping can have relevance to the appearance of super- 
conductivity in high-temperature oxide compounds. In the present paper, we investigate 
the generic occupation-dependent hopping Hamiltonians with respect to peculiarities of the 
normal state, and to the range of existence of the superconducting state. Theoretical inves- 
tigation of Cooper instability is supplemented by numeric study of pairing and diamagnetic 
currents in finite atomic clusters. We study the effect of Cooper pairing between the carriers 
and show that at certain values and magnitudes of the appropriate coupling parameters, 
the system is actually superconducting. The properties of such superconducting state are 
in fact only slightly different from the properties of conventional (low-Tc) superconductors. 
Among those we so far can only mention the change in the fluctuation conductivity above 
or near the critical temperature T^. Relaxation of the pairing parameter to equilibrium 
acquires a small real part due to the asymmetry of contraction-derived interaction between 
the quasi-particles above and below the Fermi energy. 

Oxygen atoms in the copper-oxygen layers of the cuprates (Figure 1) have simple 
quadratic lattice. We assume that pz orbitals of oxygen {z is the direction perpendicu- 
lar to the cuprate plane) are bound to the near cuprate layers whereas carriers at the Px,Py 
orbitals may hop between the oxygen ions in the plane. 

Let ti be the hopping amplitude of p^ {py) and t2 the hopping amplitude of Py {p^) oxygen 
orbitals between the nearest lattice sites in the x (y) direction in a square lattice with a lattice 
parameter a. Then the non-interacting Hamiltonian is 

Ho = -ti 4«i ~ ^2 X! ~ ^1 XI ^l^i ~ ^2 X! (1) 

<ij>x <ij>y <'i-j>y <ij>x 

where a| (a^) is the creation (annihilation) operator for p^ and correspondingly bl [hi) for py 
orbitals. The interaction Hamiltonian includes the terms 

Hi= aUj [Vniimj + W{mi + mj)] + ^ blbj [VniUj + W{ni + rij)] (2) 

<ij> <ij> 

where rij = a|aj, rrii = h\hi. This corresponds to the dependence of the hopping amplitude 
on the occupation numbers rij, rrtj of the form 

(iij) = ro(l - mi)(l - m^) + Ti [(1 - mi)mj + mi(l - m^)] + Tamimj (3) 
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and correspondingly {tij)b^^b of the same form with replaced with n^. The amplitudes 
To, Ti, T2 correspond to the transitions between the ionic configurations of oxygen: 

To: or + 0]-^0t + 0J 

n: + ^Ot + 0j (4) 
r2 : + OJ ^ or + O, 

O corresponds to the neutral oxygen ion whereas O to the single charged and O^^ to the 
double charged negative ions. Since oxygen atom has ls^2p^2s^ configuration in its ground 
state, filling of the p shell to the full occupied configuration 2p^ is the most favorable. 
Amplitudes V and W relate to the parameter tq, ti, T2 according to 

V = To-2Ti + T2, W = Ti-T2. (5) 

Assuming ti = t2 and replacing a^, hi with Oj with the pseudo-spin indices cr =|, j we write 
the Hamiltonian Eq.(l) in the form 

H = -t al„aj„ + Hu + Hv + Hw (6) 

<ij>(T 

where 



Hu 


i 


(7) 


Hv 




(8) 




<ij> cr 




Hw 




(9) 




<ij> o- 





where we also included the in-site Coulomb interaction (U) between the dissimilar orbitals 
at the same site, a can also be considered as a real spin projection of electrons at the 
site. In that case, the pairing will originate between the spin-up and spin-down orbitals, 
rather than between and Py orbitals. More complex mixed spin -and orbital- pairing 
configurations can also be possible within the same idea of orbital contraction (or expansion) 
at hole localization but are not considered in this paper. The following discussion does not 
distinguish between the real spin and the pseudo-spin pairing. The Hamiltonian, Eq.(6), 
is a model one which can not refer to the reliable values of the parameters appropriate 
to the oxide materials. The purpose of our study is rather to investigate the properties 
of superconducting transition specific to the model chosen and to find the range of the 
U, V, W values which may corresDond to superconductivity. This will be done along the 
lines of the standard BCS modetj in the weak coupling limit, U, V, W 0, and by an 
exact diagonalization of the Hamiltonian for a finite atomic cluster at large and intermediate 
coupling. 

In the momentum representation, the Hamiltonian becomes H = Hq + Hi + H2 with 

^0 = I] ^pa{,^ap^ (10) 

per 

^1 = ^ Y aii«"L/3r°/375(Pl'-P2)P3,P4)ap45ap3^ (ll) 

PlP2P3P4,0/3 7(5 
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where 



= — tcTp — /i, cTp = 2{cos p^a + cospj^a), 
and /i is the chemical potential. r°^^^ is the zero order vertex part defined as 



(12) 



U +{W + -Z/V)((Tp, + (Jp2 + (Jp3 + CTpJ r^^r^5 (5^7(5/35 - (5/37 )'^Pi+P2,P3+P4 



(13) 



where r^^ is a Pauli matrix 



(;;)■ 

For reasons which will be clear later, we separated Hy and put some part of it into the Hi 
term, while the remaining part is included in the H2 term, thus giving 



<tj>cr 



H2 = V J2 4a(^jA4a(^i^ - ^)(4<?%'? 



2' 



(14) 



with 1/ =< Hi > being the average occupation of the site. 



II. THE COOPER INSTABILITY IN THE OCCUPATION-DEPENDENT 

HOPPING HAMILTONIANS 

The Cooper instability realizes at certain temperature T = as a singularity in a 
two-particle scattering amplitude at zero total momentum. Let's introduce a function 

r(piP2, T - r') =< r^ap,|(r)a_p^(r)a_p2|(r')ap,|(r') > (15) 

where apa(r) = exp(ifr)a]3„ exp(— i/r), Opo, = exp(iJr)apQ exp(— i/r) are the imaginary 
time (r) creation and annihilation operators. At pi = — P2, Ps = — P4, the kernel of F^^^^ 
is proportional to G^g {G is one-electron Green function). We keep notation F(pp') for 
such a reduced Green function specifying only momenta p = Pi = — P2 and p' = Ps = — P4. 
By assuming temporarily \^ = 0, this Hamiltonian results in an equation for the Fourier 
transform F(p, p',^) 

F(p, p', n) = F°(p, p') - 5: F°(p, k)G^(k)G_+o(-k)F(k, p', Vl) (16) 

w k 

corresponding to summation of Feynmann graphs shown in Figure 2. In the above formulas, 
uo = {2n + 1)ttT and Q = 27rmT (n, m integers) are the discrete odd and even frequencies 
of the thermodynamic perturbation theoryl3. G'(k, cj) is a one-particle Green function in a 
Fourier representation 

G{k,uj) = -^ 



(17) 
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Diagrams of Figure 2 are singular since equal momenta of two parallel running lines bring 
together singularities of both Green functions ^(k, cj) and G'(— k, cj). 

6- vertex interaction, Eq.(8), is not generally considered in the theories of strongly- 
correlated fermionic systems. Such interaction also results in singular diagrams for p — > — p 
scattering shown in Figure 3. Since a closed loop in this figure does not carry any momentum 
to the vertex, it reduces to the average value of G which in turn is the average of the number 
operator, < a^a >. Taking into consideration of such diagrams is equivalent to replacing 
one of the nj's in Eq.(8) to its thermodynamical average u 
can be added to the renormahzed value of W, 



< al^ttia- >. Then the V term 



W 



-uV 



We will check by numeric analysis in Sec. Ill to which extent such an approximation may 
be justified. 

Solution to Eq.(16) can be received by putting 



r(p, p', O) = A{n) + Bi{n)ap + 52(0)(7p. + C(0)(7p(7p/. 
Substituting this expression into Eq.(16) and introducing the quantities 



(18) 



(19) 



we receive a system of coupled equations for A, Bi, B2,C 



f I + US0 + WSi 


USi + WS2 


^ \ 






(u\ 


WSo 


1 + WSi 





Bi 




w 








1 + USo + WSi US1 + WS2 


B2 




w 


V 





WSo 1 + WSi J 


\c ) 




[ / 



where W — W + hiyV, which are solved to give 



(20) 



U 



W^S2 „ „ W{l + WSi) 



D 



-, B1 — B2 — 



D 



where is a determinant 



D 



1 + + USx + WS2 
WSo 1 + WSi 



W^So 
D 



(21) 



(22) 



The determinant becomes zero at some temperature which means an instability in the 
two-particle scattering amplitude (F — > 00). This temperature is the superconducting tran- 
sition temperature Tf.. At Tc, Eq.(16) is singular, which means that two-particle scattering 
amplitude gets infinite. Below Tg, the finite value of F is established by including the non- 
zero thermal averages (the order parameters) 



< flpfl-p >. We first analyze the 



case of non-retarded, non-contraction interaction U ^ and after that will consider the effect 
of the occupation-dependent hopping terms, V and W . 
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A. Direct non-retarded interaction 



Neglecting contraction parameters V, W, solution to Eq.(16) reduces to 

which after the summation over the discrete frequencies reduces to the conventional BCS 
equation (at negative U) 

with nk = (exp(/5^k) + 1)^^- At finite frequency fl, Eq.(23) reduces to 

where we replaced for simplicity an integration over the Brillouin, zone / d^k, by the inte- 
gration over the energy assuming that the density of states near the Fermi energy n is flat. 
—El and E2 are the lower and upper limits of integration equal to — 4t — ^ and 4t — fi, 
respectively. Such an approximation is not very bad since most singular contribution to 
integral comes from the point = where the integrand is the largest. 

Above Tc, Eq.(25) determines the frequency of the order parameter relaxationlli0§. 
There is a small change in this frequency compared to the BCS model in which limits of 
the integration {—Ei,E2) are symmetric with respect to the Fermi energy, and small in 
comparison to ep, therefore we briefly discuss it now. 

To receive a real-time relaxation frequency, Eq.(25) needs to be analytically continued 
to a real frequency domain from the discrete imaginary frequencies iun = (2n + l)iTi 
Using the identity 

TS" - - ( f V TT ^^^'^ (26) 

where n(^) is a Fermi function n{^) = (exp(/3^) + 1)^^ gives 

tanh^ 

where 



r. = B2^exp(-^^^).ln7 = C = 0.577, (28) 

C is Euler constant. Analytic continuation is now simple: we change VL to i{uj — i5), 5 = +0, 
to receive a function which will be analytic in the upper half plane of complex u, Imuj > 0. 
The order parameter relaxation equation becomes 
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(29) 

At a; <^ Tc and T — Tc <^ T^, the real and imaginary parts of Eq.(29) are easily evaluated 
to give 

r-T.-^ + .^)A = o. (30) 



Thus, the order parameter relaxation equation at T > Tc becomes 



(l + iA)^ + rA = (31) 



where 



r = V-T^c), A^ '^V/'^ ^c (32) 

In comparison to the BCS theory in which Ei = E2 = ujd {ujd is the Debye frequency) 
and therefore A = 0, wc receive the relaxation which has a non-zero "inductive" component, 
—iXT. Typically, Ei E2 ^ sp and therefore |A| is a small quantity. It increases however 
near the low (1/ <^ 1) or near the maximal {v ~ 2) occupation where Ei or E2 become small. 
Such mode of relaxation is specific to a non-retarded (non-phonon) interaction which is not 
symmetric near ep and spans over the large volume of the k-space rather than is restricted 
to a narrow energy uj^ <^ Sp near the Fermi energy. 



B. Occupation-dependent hopping instability and relaxation 

Neglecting direct interaction, we put [/ = in Eq.(22) and receive 

-i = S^{uj) ± ^JS^{uj)S2{uj) (33) 
W 

where at finite frequency u) 

'■■^2 + tanh^ 



^ha^e (34) 



Putting (x; = we receive from Eq.(33) a transition temperature Tc- The equation has a 
solution at < 0, /i < 0, or at > 0, /X > (we assume that t > 0). Plus or minus sign is 
chosen to receive the maximal value of Tc (the second solution corresponding to smaller T, 
then, has to be disregarded since at T < Tc the order parameter will be finite and therefore 
Eqs.(20)-(22) do not apply). This gives an expression for Tc 



Tc — —\/EiE'_ 

TT • 



2exp 



'""^l 2\W\N{sp)) ^^^^ 



where n < 0, W < (second exponent is dominating the first one in the weak coupling limit 
W ^ 0). Real and imaginary parts of Sn{uj) are calculated at a; <^ Tc 
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ImSn{u)c^-^(-^) Niep) (36) 



ST, 



c 



+ ^=1 (37) 

+ V + fln^, n = 2 

Equation for A is received with a value larger than the previous one (Eq.(32)) 

Eigenvalue equation gives the p-dependence of the two particle correlator r(p, p') = 
< a|,^aLp^a-p'iap't > near 

r(P, P') =C[S2- Si{(Tp + (Tp>) + S'oCTpCTp/] . (39) 

Since C diverges at T^., this determines that order parameter becomes macroscopic at T < T^. 
Then, the pair creation operator, a^alp, will almost be a number, i.e., we may decompose 
Eq.(39) into a product 

A;Ap =< al^alp^ >< a_p/^ap/t > (40) 
and, to be consistent with the p, p' dependences, by putting ^p = ^p/ we receive 



Ap = Ci [exp{te/2)p2{0) + exp{-te/2)^So{0)j exp{tip) (41) 

where 

cos^ = -S^iO)/^JSoiO)S2iO) (42) 

and (fi is an overall phase which is irrelevant for a single superconductor but is important 
for calculating currents in multiple or weakly coupled superconductors. Therefore, system 
undergoes a pairing transition at temperature found from the Eq.(35). Since the pairs are 

charged, the state below Tc can not be non-superconducting. 

We have not calculated the Meissner response but in the following section we present 
numerical calculation of flux quantization which supports the above statement. 



III. EXACT DIAGONALIZATION OF THE OCCUPATION-DEPENDENT 
HOPPING HAMILTONIANS IN FINITE CLUSTER 

We calculate the ground state energy of a cubic system as shown in Figure 4. A magnetic 
flux $ is produced by a solenoid passing through the cube. Corners of the cube are the lattice 
sites, which can be occupied by electrons. With the inclusion of the magnetic flux, model 
Hamiltonian, Eq. 6, becomes 

H = —t al^ajcrexp{iaij) + h.c. + U^^ni^Tiii + 

<ij>a i 

+ al^aja [VniaTij^ + W{nia + nja)] e>q){iaij) + h.c. (43) 

<ij>cr 
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where 



aij = {27t/%) r A 

Jn 



dl 



(44) 



and $0 = hc/e is the magnetic flux quantum. Throughout the calculations we take t = 1. 
We start with constructing the model Hamiltonian. In a Hilbert space of one electron 



(45) 



with a basis specified as ipo = (0, 1) for the ground state (n = 0) and ipi = (1,0) for the 
excited state (n = 1). In case of states operator of annihilation a„ takes the form 





where e is the unit matrix and s is unitary matrix 





(46) 



(47) 



and stands for the Kronecker matrix multiplication. Explicitly, we have 

fli = a®e®e(8>e...(8>e 
a2 = s®a®e®e...®e 



ttN 



s . . . ® s ® a 



Thus, for example, for two states 





/O 1 






/O 1 





\ 















-1 




ai = 





1 














^0 


o) 




^0 








(48) 



These matrices, which are annihilation operators, and corresponding Hermitian conjugate 
matrices, which are the creation operators, satisfy the Fermi anti-commutation relation. 
These operators are sparse matrices with only N/2 non-zero elements, which are equal to 
±1. Next we solve the Schrodinger equation Hip = Eip. We implemented a novel algorithm 
for solving such sparse systems, which will be described elsewhere. 

The cubic cluster within the Hubbard Hamiltonian and no external flux applied to the 
system was studied previously by Callaway et. al.lii. Quantum Monte Carlo methods appli- 
cable to large systems within the Hubbard model (both attractive and repulsive), but not 
the occupation-dependent hopping Hamiltonians, are reviewed in a paper of Dagotto§. 
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A. The number parity effect 



Superconductivity reveals itself in the lowering of the ground state energy as electrons 
get paired. Therefore the energy needs to be minimal for even number of electrons n and 
will attain a larger value when n is odd. We consider a "gap" paramete r§ 

A; = E21+1 - ^ {E21 + E21+2) (49) 

as a possible "signature" of superconductivity (where Em corresponds to the ground state 
energy for m fermions). For all interaction parameters set to zero (f/ = \^ = = 0), no sign 
of pairing is observed. To check our analytic results of Sec. lib and the argument following 
Eq.(34), we calculated A above and below the half- filling (ra = 8 in case of cubic cluster). 
Below the half-filling chemical potential is negative (yU < 0) and above the half-filling it is 
positive (/i > 0). We first checked that the ^ 0+, ^ 0" and V" ^ 0+, 1/ ^ 0" 
calculation is consistent with an exact solution available for a non-interacting system of n 
electrons. 

We then test our program for the case of negative-?/ Hubbard Hamiltonian ([/ < 0, = 
0, = 0) which is known to be superconducting (e.g. Refs. 22,23). Positive-t/ Hubbard 
model does not show any sign of superconductivity, in disagreement with some statements in 
the literature@. Our calculations can not disprove the (possible) non-pairing mechanisms of 
superconductivity but these seem to be unlikely models for the problem of superconductivity 
in oxides which clearly shows pairing of electrons (holes) in the Josephson effect and in the 
Abrikosov vortices. The relation 2eV = Uuj is justified in the first caseil and flux quantum 
of a vortex is hc/2e in the secondil, both with the value of the charge equal to twice the 
electronic charge, e. 

Figure 5 shows the dependence of the ground state energy upon the number of particles 
in case of negative-?/ and positive-?/ Hubbard models assuming = and W = Such 
dependences are typical for any value of \U\. There clearly is the pairing effect when f/ < 
and there is no sign of pairing at ?/ > 0. 

Tests for pairing in the contraction V , W-modeh {V ^ U = W = Q and W 7^ 
0, f/ = y = 0, respectively) are shown in Figs. 6,7. The results are in agreement with our 
perturbative calculation of Sec. II and with its extension for the intermediate and strong 
coupling limits \V\ ^ t, \W\ ^ t. Since chemical potential is negative below the half-filling 
and positive above the half-filling, there is no pairing in the former case {W 0"^) and 
there is a sign of pairing in the latter case iW > 0), in accord with the value of the effective 
coupling constant W = W + \h'V . Similarly, for 0^ below the half-filling there is a 

sign of pairing (A 7^ 0) while above the half-filling there is no pairing. These results are 
summarized in Table 1. 

For larger values of the interaction parameters, the perturbative results do not remain 
applicable anymore. Figure 8b shows the dependence of the parity gap A on the strength 
of the interaction. From Figure 8 it is understood that the W interaction introduces a 
"signature" of pairing in a similar way as the negative-f/ interaction does. The possibility 
of the "contraction" pairing has been investigated formerly in the papersS'El. 
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B. Flux quantization 



Flux quantization is another signature of superconductivity which is a consequence of 
the Meissner effect. We also tested for the periodicity of the energy versus flux dependence 
with the period $i = hc/2e as compared to the period $o = hc/e in the non-interacting 
system0't3. Unfortunately, the even harmonics of $o-periodic dependence of the ground 
state energy (and related to it, the harmonics of the persistent current J = —dE/d^^^) 
may simulate the pairing in a non-superconductive system. Small-size (mesoscopic) system 
can mask the superconducting behavior@. Flux quantization in Hubbard Hamiltonians was 
studied formerly in Refs. 29-31. 

We first demonstrate the behavior of the ground state energy with respect to flux. Figure 
9. A characteristic feature of mesoscopic system suggests that addition of one extra particle 
to the system changes the sign of the derivative of the ground state energy with respect to 
magnetic flux at $ = 0. That is, depending on the parity of the number of particles and 
on the number of sites, system can change from paramagnetic to diamagnetic state or vice 
versa. But this behavior is not always observed for the cubic geometry studied. Except 
the sign change from n = 2 to n = 3 and from n = 7 to n = 8, no such behavior is seen. 
As mentioned above, however, the $i-periodic component of the E{^) dependence begins 
to appear at the higher value of n (Figure 9,c). For both contraction parameters equal to 
zero, i.e. W = V = 0, we observe appearance of the /ic/2e-periodic component for some 
values of U (Figure 10). Even for positive (repulsive) values of U, it is possible to see a local 
minimum appearing at $ = hc/2e (Figure 10b). This is in agreement with the authors' 
previous workslilii. But this minimum, which does not lead to an exact periodicity of the 
ground state energy with a period $o/2, should not be attributed to superconductivity, this 
is rather a characteristic behavior in mesoscopic systems. 

For U < (while W = V = 0), the expected mesoscopic behavior, that is the change of 
the sign of the slope of ground state energy at $ = 0, starts to demonstrate itself (Figure 
11). But this happens at sufficiently large absolute values of (negative) U. For other values 
of U, however, there is no such change. 

More pronounced /ic/2e-periodic components are observed with the introduction of non- 
zero interaction parameters. The role of W on the ground state energy, when both U and 
V are zero, is shown in Figure 12. Meanwhile setting both U and W to zero and observing 
the effect of the non-zero V shows that V does not play a role as significant as the other two 
interaction parameters do. There is not much difference in the behavior of the ground state 
energy upon magnetic fiux between the zero and non-zero V (for example V = —1) cases. 



IV. CONCLUSIONS 

We studied the peculiarity of electron conduction in systems in which conduction band 
is derived from the atomic shells with a small number of electrons [Ng) in an atom. Such 
materials may include oxygen [Nf. = 8) in the oxides, carbon (A^^e = 6) in borocarbides (e.g. 
LuNi2B2C) , hydrogen (A'e = 1) in some metals (e.g., Pd — H). Some materials of this 
kind are superconductors. It was argued that the Coulomb effects within the atoms strongly 
influence the inter-atom wave function overlap between the atomic sites and therefore the 
electron hopping amplitude between the sites. The phenomenology of such conduction 
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mechanism results in a novel, to the conventional solid state theory, Hamiltonians called 
the occupation-dcpcndcnt-hopping (or contraction) Hamiltonians, specified with the two 
coupling parameters W. We then attempted a study of superconductivity in such systems 
within the BCS-type approach assuming the Cooper pairing of electrons. The weak-coupling 
limit allows determination of the range of parameters V, W values and also of the in-site 
Coulomb interaction U value which show the Cooper instability. The strong-coupling limit 
was addressed by a numeric calculation on finite clusters using the novel algorithm (of non- 
Lanczos type) for eigenvalues of large sparse matrices. One of the results of this numeric 
calculation was that the positive-C/ Hubbard model, sometimes believed to be a candidate 
for high-Tc superconductivity, does not comply with the goal. 

This work was partially supported by the Scientific and Technical Research Council of 
Turkey (TUSiTAK) through the BDP program. 
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FIGURES 
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FIG. 1. Site configuration in the Cu02 plane of cuprates. Dotted fine represents the effect of 
orbital contraction/extension due to the locahzation/delocahzation of an extra hole at a specific 
site. The enlarged orbital attains the larger value of the hopping amplitude to the nearest sites. 
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FIG. 3. Cooper diagram for 6-vertex interaction, V. 
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FIG. 4. Sample configuration. The flux through the cube is produced by a solenoid. 
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FIG. 5. Dependence of the ground state energy upon the number of particles with U ^ and 
V = W = 0. (a) For U <0, the pairing effect is clearly seen, (b) For U > 0, there is no pairing 
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FIG. 6. Dependence of the ground state energy upon the number of particles with V ^ and 
U = W = 0. (a), (b) Both for V > and V < 0, around the half-fining, there is a small pairing 
effect. 
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FIG. 7. Dependence of the ground state energy upon the number of particles with W ^ and 
U = V = 0. (a), (b) Both for Vl^ > and W < 0, there is a more pronounced pairing effect below 
the half-filling. 
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FIG. 8. Dependence of the parameter A upon U for various values of W and V below the 
half-filling. 
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FIG. 10. Dependence of the ground state energy upon magnetic flux. Contraction parameters 
are both zero, i.e. W = V = 0, only the on-site interaction parameter U is nonzero. 
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FIG. 11. Dependence of the ground state energy upon magnetic flux. Comparing (a) with 
Figure 9(b) clearly shows that the change in the parity of the number of particles for the case of 
negative U values introduces a sign change in the slope of £^(<I>) at $ = 0. 
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FIG. 12. Dependence of the ground state energy upon the magnetic flux. On-site interaction 
parameter U and one of the contraction parameters, F, are zero. AU plots correspond to the 
non-zero interaction parameter W = —1. 
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TABLES 





U = W = 


U = V = 








W ^0+ 


W ^0- 


below half-filling 










(m<o) 


A = (no pairing) 


A 7^ (pairing) 


A = (no pairing) 


A 7^ (pairing) 


above half-filling 










(/x>0) 


A ^ (pairing) 


A = (no pairing) 


A ^ (pairing) 


A = (no pairing) 



TABLE I. Pairing effect for arbitrarily small values of V and W, computed by exact diag- 
onalization of the Hamiltonian. The results presented here are in complete agreement with the 
perturbative calculations. 



26 



